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ABSTRACT 

The Zeldovich approximation, order Lagrangian perturbation theory, provides a good de¬ 
scription of the clustering of matter and galaxies on large scales. The acoustic feature in the 
large-scale correlation function of galaxies imprinted by sound waves in the early Universe 
has been successfully used as a ‘standard ruler’ to constrain the expansion history of the Uni¬ 
verse. The standard ruler can be improved if a process known as density field reconstruction 
is employed. In this paper we develop the Zeldovich formalism to compute the correlation 
function of biased tracers in both real- and redshift-space using the simplest reconstruction al¬ 
gorithm with a Gaussian kernel and compare to N-body simulations. The model qualitatively 
describes the effects of reconstruction on the simulations, though its quantitative success de¬ 
pends upon how redshift-space distortions are handled in the reconstruction algorithm. 

Key words: gravitation; galaxies; haloes; galaxies: statistics; cosmological parameters; large- 
scale structure of Universe 


1 INTRODUCTION 

The large-scale structure seen in the distribution of galaxies con¬ 
tains a wealth of information about the nature and constituents of 
our Universe. Of particular interest here is the use of low-order 
statistics of this field to constrain the distance scale and growth rate 
of fluctuations, which in turn impact upon our understanding of 
dark energy and tests of General Relativity at cosmological scales 
(e.g. [Olive et aH2014) . One of the premier methods for measuring 
the distance scal^uses the baryon acoustic oscillation (BAO) ‘fea¬ 
ture’ in the 2-point function of galaxies as a calibrated, standard 
ruler (see [Olive et al.|2014| for a review). Additional information 
on the rate of growth of perturbations, which allows a key test of 
General Relativity and constraints on modified gravity (e.g. [Joyce[ 
[et al.[[20T4| and references therein), is encoded in the anisotropy 
of the 2-point function imprinted by peculiar velocities, i.e. red- 
shift space distortions (see [Hamilton| 199^ for a review). Fits to the 
distance scale using the BAO feature become significantly more ac¬ 
curate if density field ‘reconstruction’ is applied ^Eisenstein, et al.| 
[2007a[ >, but this procedure alters the signal that is used to infer the 
growth rate from redshift-space distortions. Ideally we would have 
a model which can simultaneously describe the features which are 
used to constrain distance scale and the growth of structure, since 
there is a non-trivial degeneracy between mis-estimates of distance 
and growth (e.g. Fig. 9 of [Reid et al.||201^ . A formalism which 
can be used to simultaneously describe both of these pieces of a 
redshift survey is currently not known. 


It is straightforward to form a data vector which consists of the 
correlation function pre-reconstruction on small scales and post¬ 
reconstruction on large scales. Our goal is to find a single theoret¬ 
ical framework which could simultaneously fit both parts of this 
data vectoQ Models based upon Lagrangian perturbation theory 
have been shown to do a good job of fitting the anisotropic sig¬ 
nal in the (pre-reconstruction) correlation function (see e.g. [Whitej 
[et al.[[20T5] for a recent investigation and references to the earlier 
literature). In this paper we investigate how accurately P* order La¬ 
grangian perturbation theory (“the Zeldovich approximation”) can 
be used to model the reconstructed BAO feature in the redshift- 
space correlation function of biased tracers. 

The last few years have seen a resurgence of interest in the 
Zeldovich approximation. It has been applied to understanding the 
effects of non-linear structure formation on the baryon acoustic os¬ 
cillation feature in the correlation function ( [Padmanabhan & Whit e| 
[2009[[McCullagh & Szalay[2012[[Tassev & Zaldarriaga[2012a^ and 

to understanding how “recon struction” jEisenstein, et al. 2007a 

removes those non-linearities ([Padmanabhan, White & Cohn|2009 


[Noh, White & Padmanabhan[2009[ Tassev & Zaldarriaga[2012b) . It 


has been used as the basis for an effective field theory of large-scale 
structure ( [Porto, Se nat ore & Zaldarriaga|2014^ and a new version 
of the halo model dSeliak & Vlah|2015b. It has been compared to 


“standard” perturbation theory dTassev|2014al, extended to higher 

orders in Lagrangian perturbation theory ( Matsubara|2008a|b 

jOka-j 

mura, Taruya, & Matsubara 

2011 [[Carlson, Reid & White 

20131 

Vlah, Seljak & Baldauf|2015 

1 and to higher order statistics ( Tassevj 


' And for breaking degeneracies when constraining parameters from the 
cosmic microwave background anisotropies, e.g. [Planck Collaboration! 

(Ms). 


^ Obviously, such a model would also form a good template for fitting the 
BAO peak position on its own. 
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2 White 


|2014b^ including a model for the power spectrum covariance ma¬ 
trix I Mohammed & Seljak|2014) . Despite the more than 40 years 
since it was introduced, the Zeldovich approximation still provides 
one of our most accurate models for the distribution of cosmologi¬ 
cal objects. 

The outline of the paper is as follows. Section|^contains a re¬ 
view of the salient aspects of Lagrangian perturbation theory and 
reconstruction, to fix our notation, and introduces our N-body sim¬ 
ulations. Sectionj^introduces the Zeldovich model for reconstruc¬ 
tion and compares its predictions to the simulations. We finish in 
Section]^ with an assessment of the Zeldovich approximation and 
future directions for research. 


2 BACKGROUND AND REVIEW 
2.1 Lagrangian perturbation theory 


We wish to develop an analytic description of the reconstructed cor¬ 
relation function of biased tracers in redshift space and to this end 
we use Lagrangian perturbation theoryp i BuchertjlV^ iMoutarde 


|et al.|[T9^|Hivon et al.|[T995[ [Taylor & Hamilton|1996| . In this 

section we remind the reader of some essential terminology, and 
establish our notational conventions. Our notation and formalism 
follows closely that in|Matsubara|l 2008a|b^; Carlson, Reid & White| 


pOT^ ; [^g, Reid & White| ( |2013^ ; |White| ( |2014| > to which we re¬ 
fer the reader for further details and original references. 

In the Lagrangian approach to cosmological fluid dynamics, 
one traces the trajectory of an individual fluid element through 
space and time. Every element of the fluid is uniquely labeled by 
its Lagrangian coordinate q and the displacement field 'I'fq, t) fully 
specifies the motion of the cosmological fluid. Lagrangian Pertur¬ 
bation Theory (LPT) develops a perturbative solution for *P but we 
shall deal here with the first order solution which is known as the 
Zeldovich approximation ( |Zeldovich| 1 970| ). Denote this first order 
solution as *P we have: 


We shall assume that halos, and the galaxies that inhabit them, have 
a local Lagrangian bias px(<l) = PxU[^«(q)]- |Matsubara| ( [20TT| ) pro- 
vides an extensive discussion of local and non-local Lagrangian 
bias schemes. 

This formalism makes it particularly easy to include red- 
shift space distortions. We follow the earlier papers and adopt the 
“plane-parallel” or “distant-observer” approximation, in which the 
line-of-sight direction to each object is taken to be the fixed direc¬ 
tion z- Within this approximation, including redshift-space distor¬ 
tions is achieved via 


T, ^ >p; = 


(dij + fz.zjWj 


( 2 ) 


which simply multiplies the z-component of the vector by 1 + /. 

The correlation function within the Zeldovich approximation 
then follows by elementary manipulations. Defining A = *P 2 “ ‘I’l 
and writing Fj = Ffd,) for the Fourier transform of Ffd^fq)] the 
real-space correlation function is 




{Inf 

<51 +A2d2+k-\) ^ 


' dAx dA2 

2 n 2 n 


(3) 


^ See [Bemardeau et al.] {2002} for a comprehensive (though somewhat 
dated) review of Eulerian perturbation theory. 


For convenience we define ^/.(q) = (did 2 ), f/,(q) = (diA,) = ( 62 ^ 1 ), 
and Aij(q) = (A/Aj). The vector (7,(q) = U{q)qi is the cross¬ 
correlation between the linear density field and the Lagrangian dis¬ 
placement field. The matrix A,j may be decomposed as 

= 2 [a-^ - q±(q)\ Sij + 2 [qAq) - PiK?)] (4) 

= cf-Jij + [crl-cf^\qiqj (5) 

where cr^ = |(|*Pp) is the 1-D dispersion of the displacement field, 

and 7/11 and q± are the transverse and longitudinal components of the 
Lagrangian 2-point function, 77,j(q) = ^'T,(qi)'i'j(q 2 )y In the Zel¬ 
dovich approximation these quantities are given by simple integrals 
over the linear power spectrum: 

‘A=-^2£dkPdk), (6) 

r]±{q)=:^ f dkPtik) (7) 

2n^ Jo kq 

q\\{q) ^ ^ J' 

t/(?) = -^j| dkkP2^{k)U{kq). (9) 

Taylor series expanding the bias terms and doing the Ti and 
A 2 integrations and the Fourier transform we can write 


Jo(kq) - 2 


ji{kq) 


kq 


(8) 


l+fx(r) 


/ 


d^q 


(2;r)3/2|A|‘/2 
1. 




2biU,gi + -bi^i - {b2 + biW.UjGij 

2bib2^l,U,g, + ---], 


( 10 ) 


where we have written b„ - g, = (A *),/(? - r)j and 

Gjj = (A“'),y - gigj in order to make the expressions more read¬ 
able. The generalization to redshift space follows straightforwardly 
from Eq. |^: we simply multiply hy I + f and divide the z- 
components of A“' by the same factor. 

Not all of the terms in Eq. l |10[ ( are important at the scales 
relevant for BAO. For typical values of halo bias (bi ~ 1 and 
^2 ~ 0.1), the dominant contributions to the real space correlation 
function or the monopole of the redshift space correlation function 
at r 100/i^'Mpc are from the “1”, b^^^L and -2biUigi terms. 
The other terms make up less than one per cent of the total. For 
the quadrupole of the redshift space correlation function only the 
the “1” and -2bi (/,g, terms contribute significantly (see also |Whit^ 
[20T4l Fig. 4). 


2.2 Reconstruction 


We start by reviewing the reconstruction algorithm of |Eisenstein,| 
|et al.| ( |2007a^ and its interpretation within Lagrangian perturba- 
tion thm ry jPad manabhan, White & Cohn||2009| |Noh, White &| 

|Padmana bhan|2009t . Various tests of reconstruction have been per- 

formed in |Seo et ^ ( |2010t ; |Padmanabhan et al. j2012| >; |Xu et al.| 
l |2013^ ; |Burden et aT] ( |2014^ ; Tojiero et al. ( |2014 i which also con¬ 
tain useful details on the specific implementations. 

The algorithm devised by |Eisenstein, et ar]j2007a| l is straight¬ 
forward to apply and consists of the following steps: 


• Smooth the halo or galaxy density field with a kernel .S (see 
below) to filter out small scale (high k) modes, which are difficult 
to model. Divide the amplitude of the overdensity by an estimate of 
the large-scale bias, b, to obtain a proxy for the overdensity field: 
d(x). 
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Reconstruction within the Zeldovich approximation 3 


• Compute the shift, s, from the smoothed density field in red- 
shift space using the Zeldovich approximation (this field obeys 
V • Rs = -d with the / replaced hy fib in R). The line-of-sight 
component of s is multiplied hy 1 + / to approximately account for 
redshift-space distortions. 

• Move the galaxies by s and compute the “displaced” density 
field, Sd- 

• Shift an initially spatially uniform distribution of particles by s 
to form the “shifted” density field, 6s- It is ambiguous whether this 
shift includes the factor of 1 + / in the line-of-sight direction or not. 
Including the 1 +f includes ‘linear’ redshift-space distortions in the 
reconstructed field while excluding it removed them.|Padmanabhan| 
|et al.| ( |2()12t ; |Xu et al.|P013| l and later works do not include this 
factor, but earlier papers did not distinguish between the uniform 
sample and the galaxies. We shall consider both approaches. 

• The reconstructed density field is defined as 6r = Sj - 6, with 
power spectrum Prik) oc (|dj|). 


Following |Eisenstein, et al.] ( |2007a^ we use a Gaussian smoothing 
of scale R, specifically S(k) = Throughout we shall as¬ 


sume that the fiducial cosmology, bias and / are properly known 


during reconshuction. 

Padmanabhan et al. ( 

2012|(;Xuet al.|(|2013|l; 

[Burden et al.| (|2014|l; 

Vargas-Magana et a 

. (2014J show that the 

reconstructed 2-point function is quite insensitive to the specific 


choices made, so this is a reasonable first approximation. We shall 
return to this issue in Section|4] 



k [h Mpc ^ ] 


Figure 1. The contribution to the variance of the ID Zeldovich displace¬ 
ment, per unit k,?Az- 0.5 for three different (Gaussian) smoothing scales: 

= 5, 10 and 20 Mpc (upper to lower sets of lines). For each set of lines 
the solid line is the linear theory prediction, the dashed line assumes stan¬ 
dard, 2"** order, Eulerian perturbation theory and the dot-dashed line is the 
Zeldovich approximation expanded to 2“‘* order. Except for the 5/;“* Mpc 
case all three approximations are in excellent agreement (see also Eig. 1 of 
[Padmanabhan, White & Cohn|2009y 


2.3 N-body simulations 

We use a suite of 20 N-body simulations to test how well the Zel¬ 
dovich model works. The simulations assume a ACDM cosmology 
with n„, = 0.274, Qa = 0.726, h = 0.7, n = 0.95, and erg = 0.8 and 
were run with the TreePM code described in | White |P002| |. Each 
simulation employed 1500^ equal mass (nip =: 7.6 x lO'^/t^'Mo) 
particles in a periodic cube of side length 1.5/i“'Gpc as described 
in |Reid & White) ( [20 1 1| | and |White et al.|j2011| l. Halos are found 
using the friends-of-friends method, with a linking length of 0.168 
times the mean inter-particle spacing. These are the same simula¬ 
tions and catalogs that were used in|Wang, Reid & White|(|2013|l; 
|White|j2014[ >; |White et al.| ( |2015[ > and further details can be found 
in those papers. Throughout we shall use halos with friends-of- 
friends mass in the range 12.785 < logm M;,/(/i“')Wq) < 13.085, 
with b ^ 1.7, which is one of the samples used in|Wang, Reid &| 
|White| ( |20T3l l; |White| ( (20T4l l. It has a relatively high bias, while at 
the same time a large enough spatial density to reduce shot noise to 
tolerable levels. 


The difference between (Si(q) and Spx) is higher-order in T* and so 
should be comparable to the effect of non-linearities in the densit}]^ 
Within the same approximation, solving V-Rs = -6 on the redshift- 
space field is the same as generating s(k) = -i(k/^^)(5(k).S(k) using 
the real-space field. 

To estimate the relative size of the correction to the shift terms 
coming from non-linearities in the density, we look at the con¬ 
tributions to the rms Zeldovich displacement for different (Gaus¬ 
sian) smoothing scales, R. In real space the ID displacement is 
[/ dkP(k)/(6ir^)f'\Fig. 0 shows the fractional contribution to the 
squared displacement from beyond-linear terms in P(k), computed 
from (standard) Eulerian perturbation theory or the Zeldovich ap¬ 
proximation [see Appendix|^for more details]. For smoothings of 
10 li^'Mpc or above the approximation appears to be ve ry good. We 
shall use R = 15 /z^'Mpe as our default (as used in e.g. Padmanab- 


|han et al.|2012[ [Anderson et al.|20i4{|Tojiero et al.|2014 i, unless 

otherwise specified. 

Under this approximation we compute the statistics of the dis¬ 
placed field by replacing with T* -l- s and of the shifted field by 
replacing *P with s in the formulae of 


3 ZELDOVICH RECONSTRUCTED 

With this background in hand it is now straightforward to develop 
a model for the reconstructed correlation function within the Zel¬ 
dovich approximation. 

3.1 The shift 

We will assume that the “shift” field, which is formally computed 
on the non-linear density field at the Eulerian position, x, can be 
well approximated by the negative Zeldovich displacement com¬ 
puted from the linear theory field at the Lagrangian position, q. 
This is a reasonable first approximation since such shifts are dom¬ 
inated by very long wavelength modes jEisenstein, et al.|2007b|l. 


3.2 Real space 

Let us first consider the statistics of the reconstructed field in real 
space. The reconstructed field is the sum of the displaced and 
the negative of the shifted fields of Sec. |2.2| and thus the corre¬ 
lation function has 3 terms: the auto-correlation of the displaced 
field, the auto-correlation of the negative-shifted field and the cross¬ 
correlation of the two fields: +^( w) + Each term 

will have the same functional form as Eq. Let us take each 

While the ‘shifts’ from Lagrangian to Eulerian coordinates are large in 
CDM, they are quite coherent so this approximation is not as drastic as it at 
first seems. 
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Figure 2. A comparison of the Zeldovich and N-body real-space, halo cor¬ 
relation functions pre- and post-reconstruction. The dotted line shows the 
linear theory, while the solid (dashed) line shows the Zeldovich predic¬ 
tion pre(post)-reconstruction. The squares and diamonds show the unrecon¬ 
structed and reconstructed results from the N-body simulations described in 
the text. We have used a smoothing scale of R = 15 /?“^Mpc when perform¬ 
ing reconstruction. 


in turn. The auto-correlation function of the displaced field, 
is given by Eq. 1 10 1 with —> Pi{l - S)^ when evaluating rj^ 

and rjw and one power of 1 - *S when computing U (it is unchanged 
when computing ^i). Thus for example the Ui entering the analog 
of Eq. 1 10 1 for is given by 


dkkPdk)(l-S) jdkq). (11) 


and similarly for the other terms. The auto-correlation function of 
the shifted field is similarly given by Eq. l |10^ with b[ = b 2 = 0 
(i.e. the terms in square brackets in Eq. jl0|l become 1) and Pl —» 
PtS^ when evaluating 77 ^ and 7711 which define A,j. The cross term 
between the displaced and shifted fields has P^ —» PlS{\--S) when 
evaluating 77 j^ and T/y and —> P^S when evaluating U and the 
substitutions b\ —> Ifoi, b 2 —> ^b 2 , —> 0 , —» 0 and bib 2 —> 0 

in Eq. jl0[ l, i.e. 


1+CTr) = 


= f 


chq 




'■ I ^ I O “ 1 1 / 


1. 


-jT-l) V)(r-q) |-j 


(ds)jAds)^(ds) 


( 12 ) 


A comparison of the correlation function predicted by the Zel¬ 
dovich approximation with that measured in N-body simulations 
is shown in Fig. The theory predicts that the acoustic peak (at 
r ^ 110 / 7 ^‘Mpc) is broadened by the effects of non-linear struc¬ 
ture formation and that reconstruction acts to sharpen the peak. The 
agreement with the simulations both pre- and post-reconstruction 
is quite good, as expected from the earlier work of |Noh, White| 
although in that work 2"“* order LPT was 


& Padmanabhan 


(2009 


used). While we do not have the necessary volume of simulations 
to reliably measure the peak location at sub-percent precision, we 
argue in the Appendix that the model should accurately reflect the 
manner in which reconstruction reduces the small shift in the peak 
location engendered by mode-coupling (see similar discussion in 
[Padmanabhan & White|200^. We have checked that the agreement 



s [h ^ Mpc] 


Figure 3. A comparison of the Zeldovich and N-body redshift-space, halo 
con'elation functions pre- and post-reconstruction. The dotted line shows 
the linear theory, while the solid (dashed) line shows the Zeldovich predic¬ 
tion pre(post)-reconstruction. The squares and diamonds show the unrecon¬ 
structed and reconstructed results from the N-body simulations described in 
the text. The upper set of lines are for the quadrupole while the lower set of 
lines is for the monopole except in the lower panel where the lowest dashed 
line is for the reconstructed quadrupole. Two versions of reconstmction are 
shown: (upper) with both the halos and the initially uniformly distributed 
particles shifted by the same field (lower) with the halos shifted 1 + / times 
further in the line-of-sight direction than the uniform particles. 


between the model and the simulations is qualitatively similar for 
variations in the smoothing scale between 10 to 20 / 7 “*Mpc. 


3.3 Redshift space 

Now we turn to redshift space. If we use a single field, s, to shift 
hoth the halos and the random particles (i.e. with the factor of 1 -(-/ 
in the line-of-sight direction for both) when generating s the mod¬ 
ifications to the preceeding section are small: we simply multiply 
by 1 + / and divide the z-components of A“' hy the same factor. 

The upper panel of Eig.j^shows the monopole and quadrupole 
of the correlation function in this case. The Zeldovich approxima¬ 
tion does a credible job of fitting the monopole of the redshift- 
space, halo correlation function pre-reconstruction. The agreement 
for the quadrupole moment is better than linear theory in the acous¬ 
tic peak region, but not as good as for the monopole (as expected 
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from earlier work, e.g. |White|[2()14| Fig. 2). To avoid cluttering 
the figure we have not plotted the errors on the N-body points. 
For the monopole they are generally small, but for the quadrupole 
(pre- and post-reconstruction) they are significant. In the acoustic 
peak region the typical error on is 3 - 5 /i““Mpc^ and the er¬ 
rors are highly correlated. Post-reconstruction the results for both 
multipoles of the correlation function are qualitatively similar: the 
reconstructed multipoles are closer to the linear theory than the 
evolved ones and the agreement with the N-body simulations in 
the region of the acoustic peak (j =: llOA^'Mpc) is quite good. 
Unfortunately the errors on the quadrupole from the N-body sim¬ 
ulations are too large to see whether the predicted shift from the 
pre- to post-reconstruction shape near the acoustic peak is home 
out in simulations. If pushed to smaller scales the model starts 
to depart significantly from the simulation results, no doubt be¬ 
cause the Zeldovich approximation does not accurately capture the 
anisotropies in the displacement/velocity field on smaller scales 
(see |White|2014| for further discussion). There is weak evidence 
that the Zeldovich approximation agrees better with the N-body 
simulations for the quadrupole moment after reconstruction than it 
does before. Increasing the smoothing scale (to 30/i“*Mpc) leads to 
similar agreement between the simulation and model, but reduces 
the sharpening of the peak by reconstruction. Reducing the smooth¬ 
ing scale to lO/j^'Mpc gives results very similar to those shown in 
Fig.|3] 

An alternative formulation does not include the factor of 1 -t / 
in the line-of-sight shift for the initially uniformly distributed par¬ 
ticles. This acts to reduce the effects of redshift-space distortions 
in the reconstructed density field. In this case the factors of I -t / 
are omitted entirely when computing the shift-shift auto-correlation 
function, and only one power of 1 -t / is included in and no 
factors of 1 -)- / in U in the cross-correlation of the displaced and 
shifted particles but the rest of the terms remain unchanged. This 
is shown in the lower panel of Fig. and the level of agreement 
between the theory and the simulations is similar to that in the up¬ 
per panel. Note in the lower panel the quadrupole is significantly 
reduced in both the model and the simulations, indicating that we 
have removed most of the effects of linear redshift-space distor¬ 
tions, but it is not reduced entirely to zero (earlier investigations of 
reconstruction in simulations either did not include redshift-space 
distortions or presented only the monopole statistics). Again the 
numerical errors from the N-body simulation are not negligible, 
but the overall trends are clear. The agreement between the simula¬ 
tions and the model in the monopole is no longer as good on scales 
smaller than the acoustic peak as it was in the upper panel. 

Comparing the upper and lower panels of Fig. [^suggests that 
the errors in how the Zeldovich approximation models reconstruc¬ 
tion partially cancel if both the galaxies and initially uniformly dis¬ 
tributed sample of particles are shifted by the same field. In this 
case the agreement between the model and simulations in both 
the monopole and quadrupole moments of the correlation func¬ 
tion above 90 /i^'Mpc is quite encouraging. If only the galaxies are 
shifted by an additional factor of 1 -(- / in the line-of-sight direction 
the reduction in the quadrupole moment is qualitatively reproduced 
by the model but the well-known inaccuracies in the halo veloc¬ 
ity field cause a significant over-estimate of the monopole even 
at 90/i^'Mpc. If the Zeldovich approximation is to be used as a 
template for fitting the reconstructed BAO feature, it would be bet¬ 
ter to implement reconstruction on the data using the ‘both shift’ 
formalism. If the behavior of the model is improved because the 
‘both shift’ formulation reduces sensitivity to small scales (where 
the model does less well) then this formulation may be less sensi¬ 


tive to small scales in the data as well and potentially more robust. 
Such an investigation is outside the scope of this work. 


4 DISCUSSION 

The goal of this paper was to investigate a model for the re¬ 
constructed, redshift-space correlation function of biased tracers 
within the framework of Lagrangian perturbation theory. In prin¬ 
ciple such a model can be combined with other models within 
the same framework to fit a combination of data such as recon¬ 
structed BAO and redshift-space distortions, for example by fitting 
a data vector which consists of pre-reconstruction multipoles below 
j ~ 90/i“'Mpc and reconstructed multipoles above i- 90/i^'Mpc. 


Previous work {[Padmanabhan, White & Cohn||2009[ 

Noh, 

White & Padmanabhan 2009|l developed the iPT formalism of 

Mat- 

subara ( 2008a|b i to reconstruction in real space and made compar- 


ison to N-body simulations. In this work we have specialized to 
lowest order in LPT, i.e. the Zeldovich approximation, but avoided 
some of the perturbative expansions inherent in iPT, extended the 
model to include redshift-space distortions and compared to a 
larger set of N-body simulations. 

The Zeldovich model performs very well, in comparison to N- 
body simulations, for the real-space correlation function of halos 
both pre- and post-reconstruction. In redshift space the monopole 
moment of the correlation function is well reproduced, and the 
quadrupole moment is consistent near the acoustic peak. Post¬ 
reconstruction the model correctly reproduces the sharpening of 
the acoustic peak and the modification of the quadrupole, but the 
quantitative agreement is not as good as in real space. The range 
of scales over which the model and the simulations agree depends 
upon how the reconstruction algorithm is implemented, with best 
agreement if both the ‘displaced’ and ‘shifted’ fields are shifted by 
the same amount. 

We have concentrated on developing and validating the Zel¬ 
dovich approximation for reconstruction, assuming that implemen¬ 
tation details, survey non-idealities and misestimates of the various 
parameters in reconstruction introduce effects that are subdominant 
to the statistical errors. This is likely true for the current generation 
of surveys (e.g. [Padmanabhan et ah |[2012[ [Anderson et al.||2014| l 
but may need to be revised for future surveys. One possibility is to 
rerun reconstruction, and recompute the 2-point statistics, for each 
cosmology whose likelihood is being evaluated (in which case the 
fiducial cosmology, bias and growth factor will be self-consistently 
included). This is extremely expensive, computationally. For small 
variations in parameters it may be possible to develop a linear re¬ 
sponse model for the 2-point function, or an emulator. Alterna¬ 
tively, an obvious direction for development is to model mises¬ 
timates of b, f and the fiducial cosmology within the Zeldovich 
approximation. This adds significant complexity to the calculation 
and obscures the main points of this paper, but may be a more com¬ 
putationally efficient method of proceeding when fitting data. As a 
side benefit it could allow an analytic understanding of the manner 
in which such assumptions impact the inferences. We defer such 
development to future work. 

I would like to thank Shirley Ho for helpful comments on an 
earlier draft. This work made extensive use of the NASA Astro¬ 
physics Data System and of the astro-ph preprint archive at 
arXiv. org. The analysis made use of the computing resources 
of the National Energy Research Scientific Computing Center. 
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APPENDIX A: ZELDOVICH VS. EULERIAN PT 


Here we briefly discuss the second order contributions to P(k) in 
(standard) Eulerian perturbation theory and in the Zeldovich ap¬ 
proximation. In the latter case it is possible to write down an ex¬ 
pression for P(k) to infinite order, but here we shall focus on the 
2"‘* order contributions. 

In both cases the second order contribution is the sum of 

/ t3 

^F2(p, k - p)F2(-p, p - k) Pdp)PL(\^ - Pi) 

(Al) 

and 


P‘^^-^\k) = 6Pdk) J £^P3(k,p,-p)PL(p) (A2) 

where F„ are the well-known perturbation theory kernels (e.g.|Go-| 
|roff et al.|1986[[Bernardeau et al.|2002t . 

Eor the Zeldovich approximation we have ^Grinstein & Wise| 

[TWtI i 



where k = p,. Thus 


F2(p, k - p)F2(-p, p - k) 


1 

4 


k • p(k^ - k • p) 

p2|k-pp 

P^(P - rf 


4(1 - 2yur -I- 


(A3) 


(A4) 

(A5) 


where we have written p -k = and p = kr. For standard perturba¬ 
tion theory ([Bernardeau et al.|2002^ 


, 5 1 Pi P2 

F2(pi,P2) - = + - 

7 2 pip2 


(Pi_ 2 (pi ■ p2)^ 

\P2^pJ^^ p]pI 


thus 


Flip, k - p)F2(-p, p - k) 


1 {Ip + 'ir-mp-rf 

T% F{\-2pr + Ff 


(A6) 


(A7) 


In both cases the integral over the azimuthal angle is trivial, and we 
are left with the p and r integrals: 


P'^-^>(/S:) = — I drP(kr) dpP(k3/l+F-- 2rp') F^(r,p) 

(A8) 

where we have written F^ipp) as a short-hand for the expressions 
in Eqs. l |A5|A7[ . 
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Reconstruction within the Zeldovich approximation 1 


For we need to evaluate p, -p). In the Zeldovich 
approximation we have 

F3(k,p,-p) = -;i^ (A9) 

3! F 

while for standard perturbation theory the expression involving the 
symmetrized form of F 3 is quite lengthy and won’t be reproduced 
here. Performing the azimuthal integral we then obtain the well 
known result for pd,3) (jjg Zeldovich approximation: 

^ -/dPdk)J^ 0 PPp) (AlO) 

while for standard perturbation theory 




Pdk) 


^ r 

- Jo 


drPdkr) 


IOOStt- 

'iid - Ifild + 2 ) 


^ - 158 + lOOr^ - 42r^ 
r- 


In 


1 +r 


1 -r 


(All) 


It is well established that Lagrangian perturbation theory, and 
the Zeldovich approximation, accurately describe the broadening 
of the acoustic peak. At this point it is also straightforward to un¬ 
derstand the origin of “shifts” in the BAO peak position due to non¬ 
linear evolution (see also |Crocce & Scoccimarro|2008|[Padmanab-| 
|han & White|200^ for discussion). Writing the convolution term 
in 

^F2(p,k-p).5z.(p)<Sz,(k-p) + --- (A12) 

in configuration space we have for standard perturbation theory 
( [Bouchet et al.|1992[|Sherwin & Zaldarriaga|2012^ 

S^6i. + ^5l + '¥-W6i, + ^T^ + --- (A13) 

where T represent (traceless) shear terms and the T • Vd term (from 
the Pi • P 2 term in Eq. |A6^ is largely responsible for the shift of 
the peak. In the Zeldovich approximation the expansion to second 
order is 

S = 6l+^61 + 'V + ■■■ . (A14) 

Note that the shift term is the same, but the growth and 
shear/anisotropy terms are slightly different (these terms match if 
we include the 2"^' order Lagrangian kernel, i.e. use 2LPT rather 
than Zeldovich). This suggests that the Zeldovich approximation 
should approximately predict the small shift in the acoustic peak 
due to mode-coupling as structure goes non-linear and the dimu¬ 
nition of this effect due to reconstruction. Further discussion and 
comparison of Eulerian and Lagrangian theories in the special case 
of one spatial dimension can be found in|McQuinn & White|(|2015l>. 
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